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DEFINITION OF SYMBOLS 


Symbol 

a 



dm 



0 

a 


n = 

= (y P >v 


Definition 

relative acceleration of a mass element as observed in moving reference frame, 
coordinate transformation from the main body system to the appendage system, 
specially defined Coriolis matrix. Format is shown in body of report, 
element of mass. 

angular momentum of rotating appendage about its mass center, 
angular momentum of vehicle, including non-rotating appendages, 
angular momentum of the total system about its mass center, 
inertia matrix for the total system, a time varying quantity, 
inertia matrix for appendage 

position vector of appendage mass center relative to the composite mass center, 
position vector of appendage mass element relative to the appendage mass center, 
vector sum of all external torques about the composite mass center, 
position vector of mass element relative to composite mass center, 
position vector of appendage mass element relative to composite mass center, 
absolute acceleration of mass element. 

absolute acceleration of the origin of the moving reference frame, 
relative velocity of mass element observed in the moving reference frame, 
rotation angle of appendage based on orbital position, 
rotation angle of appendage based on sun position. 

angular velocity vector of moving reference frame (main body) with respect to inertial 
space. 

angular velocity vector of appendage reference frame with respect to the moving reference 
frame. 

angular acceleration of the main body relative to inertial space, 
angular acceleration of an appendage relative to the main body. 


IV 



(ft)j 

(H) v 


E 

l 

& 

C 


time rate of change of vector H relative to inertial space, 
time rate of change of vector H relative to the main body, 
unit dyadic. 

inertia dyadic for total system, 
inertia dyadic for appendage. 

Coriolis dyadic. 


v 



TECHNICAL PAPER 


SPACE STATION ROTATIONAL EQUATIONS OF MOTION 
I. INTRODUCTION 


This report derives the rotational equations of motion for a large space structure having append- 
ages that can rotate relative to the main body. In particular, the equations were formulated for use in 
the Space Station Attitude Control and Stabilization Test Bed. The solar arrays and thermal radiators 
of the Space Station are required to maintain a specific alignment with the sun, whereas the main body 
of the Space Station rotates at orbital rate in an Earth pointing attitude. Thus the angular velocity of 
the appendages relative to the Space Station is approximately opposite to the orbital rate. 

Two methods are given for the derivation of the rotational equations of motion. The first 
method uses D’Alembert’s principle whereas the second method is based on the angular momentum con- 
cept. It is shown that the two formulations are dynamically equivalent although the resulting differential 
equations are significantly different concerning their mathematical structure. It should be noted that 
throughout the report the terms appendage, solar arrays or panels and thermal radiator are used 

interchangeably. 


II. DERIVATION OF ROTATIONAL EQUATIONS OF MOTION 


A. Method 1 : D'Alembert's Principle 

In this method, the rotational equations of motion are derived by applying D’Alembert’s principle. 
The equations of motion are presented in a form suitable for attitude and stabilization studies of the 
Space Station. Its Power Tower configuration is schematically shown in Figure 1. The rotational equa- 
tions will be written in terms of a moving reference frame, whose origin is located at the mass center 
of the total system and which rotates with the main body. In Figure 2 the XYZ system is fixed in 
inertial space and the xyz system is the moving reference (body) frame which translates and rotates 
relative to it. Notice that this body frame is, in general, not fixed in the main body but has a transla- 
tional motion relative to it because of the change in location of the composite mass center relative to the 
main body due to internal moving components. However, the rotational rate of the body frame is 
identical to that of the main body at all times. 

Now, D’Alembert’s principle states that the vector sum of the moments of all inertial and external 
forces relative to any point is zero. Taking moments about the composite mass center of the Space 
Station, the following vector equation is obtained: 


L + f r x (-R dm) = 0 . ^ 

The absolute acceleration R of each mass particle dm can be expressed in terms of quantities measured 
in the rotating reference frame by the well-known relation: 




z 



Figure 2. Inertial and moving reference frames. 
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R = R Q + (0 x r) + x (fl x r) + 2 (fl x v) + a 


( 2 ) 


The only components moving relative to the main body considered here are the solar panels. Their 
relative motion consists of angular rotations relative to the main body. The solar panels are assumed 
to be perfectly symmetric with respect to their axes of rotation. As a consequence, the mass center of 
the solar panels always remains at rest relative to the main body during their rotational motion. This 
restriction can be easily lifted if a higher fidelity model should be required later. The relative velocity 
v and relative acceleration a of a solar panel mass element are then given by 


Y = y p x J p . w 

a = cbp x l p + w p x (qjp x i p ) . ( 4 ) 

The position vector jp of a solar panel mass element with respect to the composite mass center can be 
expressed as: 


2p -o + -p 


(5) 


where J 0 is the position vector of the solar panel mass center relative to the composite mass center. 

Using equations (3) and (4), the five terms of equation (2) are substituted into equation (1). 
Now each term will be individually evaluated. 


TERM I: 


y*r x R Q dm = J’rdm x R Q = 0 


(I) 


The vanishing of this term is caused by the fact that the origin of the body frame was chosen to coincide 
with the composite mass center at all times, or symbolically 


f 


r dm = 0 


( 6 ) 


TERM II: 


J*[l x (fi x r)] dm = / • (2 . 

This term is brought into a vector-dyadic form by introducing the moment of inertia dyadic 
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I = ~ (r r)] dm 


where E is the unit dyadic and the integration extends over the total system mass. Thus Term II is 
obtained. 

TERM III: 


Jr x [J2 x (S2 x r)] dm = n x J [ r x (£2 x r)] dm = S2 x / • J2 


To get Term III, use is made of the vector identity: 


r x [fi x (fl x r)] = SI x [r x @ x r)] 


Term III is sometimes referred to as Euler coupling. But since it is the inertial torque obtained by sum- 
ming up the moments of the centrifugal forces acting on each mass element it is more appropriate from a 
dynamic point of view to call it the centrifugal torque. Notice that the integration of Terms II and III 
extends over the total system including the articulated appendages. Moreover the moment of inertia 
dyadic relates to the instantaneous system configuration, and thus, changes with time whenever the 
appendages are in motion relative to the main body. 

The last two terms of equation (2) reflect the influence of the rotating appendages on the attitude 
dynamics of the Space Station. When inserted into equation (1) the corresponding integration extends, 
of course, only over the appendage mass. As a reminder mathematical quantities pertaining to the 
appendages are identified by the subscript p. 

TERM IV: 


2 |r n x [fl x (cop x ip)] dnip = 2 Q. ■ C n x gj 


=P -P 


As a preliminary step, equation (5) is used to obtain the result that 


ip dm = 0 . 


Using equations (3), (4), (5) and (9), Term IV becomes 


2 J't x (J2 x v) dm p = 2 Jr p x [O x (qj p x i p )]dm p = 2 J x i p ) x [£2 x (cp p x i p )] dm p 


= 2 l ^ x [£2 x (cp p x i p )]dm p . 
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Applying the vector triple product rule 


n x (<y p x jp) = (S2 • i p ) <±> p - (0 • y p ) ip 


(ID 


equation (10) can be reformulated to obtain 

2 fl p x m • i p ) y p - (£ • Up) ipldm p = 2^(5 • l p ) (ip x qj p ) dm p . (12) 

Finally the last term of equation (12) can be converted to vector-dyadic form by introducing the Coriolis 
dyadic 


Qp = /a P ip ) d m P = (0 - 5 * TR i P )i-i P • (13) 

where TR is the trace operator. 

Note that Term IV represents the Coriolis torque (i.e., the moment of the Coriolis forces) acting 
on the appendage about its mass center. 

TERM V: 


f (Ip x a) dm p =/ p • y p + oj p x 7 p • 


(V) 


Using equation (5) and the vector identity (8), Term V is obtained. 

This last term can also be converted to vector-dyadic form by defining the moment of inertia 
ip of the appendage about its center of mass 
using equation (4), the following equation is obtained: 


dyadic / p of the appendage about its center of mass. Taking steps analogous to the previous ones and 


f(£p x a) dm p = J r p x [(oj p x l p ) + w p X (to p x Jp)] dm p 

-/a o+Jp) x [(co p x I p ) + y p x (oj p x I p ) dm p . 


(14) 


Combining Terms (I) to (V) yields the rotational equation of motion in vector-dyadic notation as: 


/ • n + J2 x / • n + 2 £ • £ p x co p +i p • cbp + y p x i p • y p = L . 


(15) 
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Notice that the first two terms represent the well known Euler equations for rigid body rotation, 
whereas the last three terms reflect the contribution of an appendage rotating relative to the main body 
with angular velocity top. If there are N rotating appendages, these terms are simply replaced by their 
appropriate sums, i.e. 


N N N 

J • $2 + x / • £2 + 2 ^ Q • C p x q; p ip • cjp + ^ ojp x /p • y p = L . (16) 

p=l p=l p=l 


B. Method 2: Angular Momentum Principle 

The angular momentum principle states that the inertial time rate of change of the angular 
momentum of a dynamic system about its mass center is equal to the sum of all external torques 
(moments of the external forces) about its mass center. In equation form, this is given by: 


L = (H t )j 


(17) 


The letter subscript outside the parenthesis indicates the reference frame in which the time rate 
of change of the vector Hj is observed. Using the known identity for vector differentiation the rotational 

equation of motion (17) can be written in the form 


L — (H'j')y + fi x Hj 


(18) 


It will now be shown that this equation is dynamically equivalent to the rotational equation of 
motion (15) obtained by applying D’Alembert’s principle. 

Before continuing with the equation development, a comment is in order regarding the usage of 
equation (17) to describe the rotational dynamics of a system. This form of equation has been used 
extensively in the past at MSFC for rigid body computer simulations in which Control Moment Gyros 
(CMG’s) are used as means of attitude control. The angular momentum of the total system is then 
defined as the sum of the angular momentum of the main body and the CMG’s about the mass center of 
the total system. Integration is performed directly on the term (H-p) v to obtain the total angular 

momentum. The desired angular velocity of the main body can then be calculated by subtracting the 
CMG angular momentum from the system angular momentum and premultiplying the resulting main body 
angular momentum by the inverse moment of inertia matrix. Because the rotating appendages of the 
Space Station are dynamically equivalent to CMG’s the simple form of equation (18) is directly applicable 
here. The present method allows, therefore, to use the existing software with only minor modifications. 
Continuing with the objective to show dynamical equivalency of the two methods, the total angular 
momentum is defined as: 


HT = H v + H p =1 ■ f2 + 7 p • top 


(19) 
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Differentiation of equation (19) with respect to time yields 

(H t ) v = 7 • (£2) v +/ p • (co p ) v + (b v • a + (7p) v • y p . (20) 

Since the rotating appendages are the only components which cause the moment of inertia of the system 
to change with time: 

(/) v = (/ p ) v • (21) 

Regarding the differentiation of the two angular velocities £2 and y in equation (20), it is noted 

that 

(£2) y = (£2)j and (y p ) v = (y p ) p . (22) 

The time rate of change of the moment of inertia dyadics can be expanded in a similar way as 
was done for the vector differentiation in equation (18) by using the known identity for dyadic differ- 
entiation in a rotating reference frame (Appendix B). 

(I>A = 0B + y BA x l ~l x ^BA • (23) 

where yg^ * s t ^ ie an gular velocity of frame B relative to frame A. Using equation (23), write: 

(ip\ = y P x i P -ip x y P (24) 

+ 

since (Z p ) p = 0. Substituting equations (19), (20), (21), and (22) into the rotational equation (18) 
we obtain: 

L = / • (£2) v + / p • (y p ) p + (i p \ • (£2 + <±>p) + £2 x (7 • £2 + 7 p • y p ) . (25) 

Then substituting equation (24) into equation (25) yields 

L = / • (6) v +i p • («p) p + (tf p x 7 p - I p x co p ) • (£2 + y p ) + £2 x (7 • £2 + 7 p • y p ) . (26) 

Since 7 p x y p • y p = 0, the equation simplies further to the final form: 
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( 27 ) 


l = / • (n> v + / p • (aj p )p + n x / ■ n + WpX/p - y p + [w p x i p • n - g p x u p ) • n 

+ £2 X / • COp] 


For dynamical equivalency of equation (27) and equation (15), the three terms in the square 
brackets of equation (27) must be identical to the term 2£2 • x (Op. The proof of this identity is 

given in Appendix A. Thus, it has been demonstrated that the two methods of deriving the rotational 
equations of motion are equivalent in a dynamic sense. 


III. VECTOR DYADIC TO MATRIX TRANSFORMATION 


For the purpose of computer programming, all vector dyadic equations have to be transformed 
to matrix equations. One of the major advantages of the vector-dyadic notation is, of course, its 
independence from any particular coordinate system. By contrast the matrix formulation requires the 
definition of a concomitant coordinate system. 

The conversion from vector-dyadic to matrix notation is based upon the isomorphism (i.c., one- 
to-one correspondence) which exists between the two. In order to obtain these correspondences , two 
vectors are considered in an orthogonal coordinate system with unit vectors e 1# e 2 , e 3 along the axes 
such that 


a = aj e^ + a 2 e 2 + a 3 e 3 , b - bj ej + b 2 e 2 + b 3 e 3 


(28) 


where aj and bj (i = 1,2,3) are the components of the two vectors along the three axes. The easily veri- 
fiable isomorphisms are obtained: 


Scalar Product: 


a • b o a 1 b . 


Vector Product: 


a x b o a" b 


where 


a = 


0 -a 3 

a 3 0 

-a 2 aj 


a 2 

" a l 

0 
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The inertia dyadic is similarly expressed in terms of the specified orthogonal coordinates as: 


/= x n + Ij2 5l 52 + ! 13 5l £3 + T 21 -2 ®1 + r 22 52 52 + *23 §2 53 + ! 31 S 3 Si 

+ I 3 2 e 3 e 2 + I 33 e 3 e 3 . (30) 

With this the following isomorphisms can also be established: 

Dyadic Product (A): 


/• w I w 

where I is the moment of inertia matrix and u is a column matrix having the vector components of 
co as elements. 

Dyadic Product (B) : 

03 • / o op I 

where co^ is the transpose column matrix. 


IV. APPENDAGE ORIENTATION KINEMATICS 


It was mentioned already that the matrix form of the equations of motion require the adoption 
of a specific reference frame (coordinate system). That is to say, all matrices and vectors appearing in 
the various terms of the following matrix equation have to refer to one and the same coordinate system, 
which in our case will naturally be the main body reference frame. Suppose that the direction cosine 
matrix specifying the relative orientation of an appendage to the main body is called Ap. Then the 

transformation of a column matrix (vector) x in the main body frame can be transformed to a corres- 
ponding vector x' in the appendage frame by the matrix relation 

x' = A p x . (31) 


Likewise the transformation of a matrix I is accomplished by 


I' = A 


P 



(32) 
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Here the primed quantities refer to the appendage coordinate system and the unprimed quantities refer 
to the main body coordinate system. 

For the following matrix equation, however, the required transformation is just reversed, namely, 
the vectors and matrices pertaining to the appendages have to be transformed to the main body frame. 
The following relation is used for the vectors: 

x = A p T x' (33) 


and the following relation is used for the matrices: 


I ApT I Ap 


(34) 


The direction cosine matrix A p is defined in this report by the following sequence of rotations 
(Fig. 1). 

1) A clockwise rotation a about the solar panel y-axis. The associated transformation matrix is 


cos a 


A y (a) = 


0 

sin a 


0 

1 

0 


-sin a 
0 

cos a 


(35) 


2) A clockwise rotation (3 about the solar panel x-axis. The associated transformation matrix is 


A X W) 


1 0 0 

0 cos P sin P 

0 -sin P cos p 


(36) 


The total direction cosine matrix is then 
A p = A x (0)A y (a) - 


cos a 

0 

- sin a 

sin a sin P 

COS P 

cos a sin j3 

sin ex cos P 

-sin P 

cos a: cos /3 


(37) 


( 38 ) 
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It is also necessary to express the angular velocities a and /3 of the solar panels in terms of com- 
ponents along the three axes of the solar panel coordinate system. This can be done by adding these 
angular velocities vectorially using the above direction cosine matrices. The result is given in matrix 
form as 


(° 

cu p = A x (/3) A y (a) I a 
' 0 


+ A x (j3) 



(39) 


From this equation the angular velocity components along the three axes of the solar panel coordinate 
system are obtained as 

cu x = )3 , cjy = a cos(/3) , co z = -a sin(/3) . (40) 


Having established these kinematical relationships, we can now transform the vector-dyadic equation (27) 
into a matrix equation using the isomophisms described in Section III as: 


i S + 3 1 n - £a p t S p c p * a p n + £ A p T i p y p + £ a p T “p ‘ p = t + L a p T ip 

P P P P 

(41) 

where L is the total torque acting on the main body and Lp is the total external torque acting on the 

p-th solar panel, properly transformed to the main body reference frame. The matrice C p * is defined 
as (cf., equation 13): 


> + ^zp J xp 

-2 J 
^ J xyp 

-2 J 
* xzp 

-2 J 
* xyp 

J xp + J zp - J yp 

-2 J 
z yzp 

-2 J 
* xzp 

-2 J 

yzp 

J xp + ^yp ” ' 


(42) 


where J vri , J,,^, J 7r . are the moments of inertia of the p-th solar panel about the mass center and J v „ n 

yy 

J XZ p, Jy Z p are the corresponding products of inertia. 


V. CONCLUSION 


It was shown that the two methods for deriving the rotational equations of the Space Station 
with articulated appendages are dynamically equivalent. It is rather surprising that the resulting set of 
ordinary first Order differential equations having such vastly different mathematical structure does indeed 
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accurately represent the very same dynamical system. The equations based upon the angular momentum 
method are quite simple in structure and many existing rigid body computer programs can easily be 
modified to include the rotating appendages. On the other hand, the equations based upon D’Alembert’s 
principle provide more insight into the dynamic origin of their individual terms, which can be of con- 
siderable benefit in trying to understand the dynamic behavior of the system. Moreover, further refine- 
ments such as mass center offset of the rotating appendages and their flexibility require only additional 
terms and no major restructuring of these equations. Ultimately, however, the choice of selecting the 
proper set of dynamic equations is largely a matter of personal predilection and experience with existing 
computer codes. 
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APPENDIX A 


The purpose of this appendix is to prove the identity 

2 S2 • Cp x Up = w p x 7 p • n - 7 p x y p • x 7p • y p . (A-l) 

As a preliminary step, it must first be shown that 

/p x Up • = /p • Wp x . (A-2) 

To this end, it is assumed for a moment that the inertia dyadic is defined by the dyadic product of two 
vectors a and b as 

7p = a b . (A-3) 

Substituting in equation (A-2) yields 


(a b) x yp • = a (b x Wp • Q) . 

(A-4) 

Since the last term in the bracket represents a scalar triple product, the 
operation can be interchanged, such that 

cross operation and dot 

a (b x Wp • f2) = a (b • ojp x £}) 


= a b • Up x Q, . Q.E.D. 

(A-5) 

To prove identity (A-l) it is best to return to the original definition of the moment of inertia. 
The first term on the right side of equation (A-l) can then be written as 

ojp x 7p • = J Wp x [r x (fl x r)] dm 


= J y p x [(r • r) n - (r • S2) r] dm 


= f (r • r) (Up x J2) - (r • tt) (Up x r) dm . 

(A-6) 


The last term of the right hand side of equation (A-l) is transformed likewise as: 
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ft x ip ' = f (r • r) (fi x ojp) - (r • w p ) (ft x r) dm . 

The second term of the right hand side of equation (A-l) is transformed as follows: 

- 7 p • (<9 p x ft) = / p • (Q x y p ) = J ^r x [(ft x w p ) x r] j dm 

= -J jr x [r x (ft x (o p )] j dm 
= - J*|r x t(r • co p ) ft - (r • ft) <y p ] J dm 

= ■ 0) (£ x y p ) - (i • y p ) (L x £?)]} dm 


(A-8) 


Adding the bottom lines of equations (A-6), (A-7), and (A-8) results in the cancellation of several terms 
with the final result: 


y p x /p*$Wp x 


ft + ft x 7 p • Wp = 2 J (ft • r) (r x w p ) 


dm 


(A-9) 


The Coriolis term of the left hand side of equation (A-l) is now transformed, going back to the defini- 
tion given in Term IV. Then 


2 ft • C p xy p = 2y*rx [ft x (co p x r)] dm 

= 2 J r x [(ft ■ r) y p - (ft • y p ) r dm 
= 2 J (Q ■ r) (r_ x cop) dm . 


(A- 10) 


It is seen that equation (A-10) is identical to equation (A-9). This completes the proof. 
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APPENDIX B 


This appendix develops the identity for the dyadic differentiation in rotating coordinate systems. 
Suppose that the dyadic / is viewed by an observer in a coordinate system A and also by another 
observer in a coordinate system B, which rotates relative to system A. Let the rotation rate of system B 
relative to system A be defined by the angular velocity vector y AB . Then the following identity holds 
for the time derivative of a generic dyadic /: 


0A “ ^B + ^AB x i~l x ^AB 


(B-l) 


To show this to be true let the dyadic / be defined by the dyadic product of the two vectors a and b, 
i.e., 


/ = a b 


(B-2) 


Then the time derivative of / as viewed in system A is: 

(t>A = 0A b + a (b) A . (B-3) 

Using the known identity for the vector differentiation in rotating coordinate systems 

(a) A = (a) B + <y AB x a (B-4) 

equation (B-3) can be written' as: 

0A = l(S>B + <°AB x al b + a [(b) B + oj ab x b] 

= [(a) B b + a (b)g] + y AB xab-ab x y A g . (B-5) 

Note that the square bracket term of (B-5) is nothing other than the time derivative of the dyadic = a b 
as viewed in system B. Thus the final result is obtained: 

0A = 0B + ^AB x l~l x ^AB ^ B ' 6 ^ 

This completes the proof. 
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